Dimensional-Crossover-Driven Mott Transition in the Frustrated Hubbard Model 
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We study the Mott transition in a frustrated Hubbard model with next-nearest neighbor hopping 
at half-filhng. The interplay between interaction, dimensionality and geometric frustration closes 
the one-dimensional Mott gap and gives rise to a metallic phase with Fermi surface pockets. We 
argue that they emerge as a consequence of remnant one-dimensional Umklapp scattering at the 
momenta with vanishing interchain hopping matrix elements. In this pseudogap phase, enhanced 
d-wave pairing correlations are driven by antiferromagnetic fluctuations. Within the adopted cluster 
dynamical mean-field theory on the 8x2 cluster and down to our lowest temperatures the transition 
from one to two dimensions is continuous. 
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The relative importance of spatial versus local fluctua- 
tions in the understanding of the Mott transition [1 can 
be tuned with dimensionality. Starting from the high 
dimensional limit, experimental studies on V2O3 indi- 
cate that exactly as a usual gas-liquid transition, the 
three dimensional bandwidth-controlled Mott transition 
belongs to the conventional Ising universality class: it 
is a first-order transition below the critical endpoint at 
Tc 450 K and affects solely the charge sector [2 . In 
contrast, Ising universality class has been ruled out in 
two-dimensional (2D) organic salts /^-(BEDT-TTF)2X 
(BEDT-TTF: bis(ethylenedithio)tetrathiafulvalene, X: 
monovalent anion) with a much lower critical point Tc — 
40 K [31, '4|. In this case geometric frustration, inherent 
to the triangular lattice, strengthens spin fluctuations 
which in turn affect the nature of the transition. The 
unconventional character of the quantum criticality in n- 
(BEDT-TTF)2X has been confirmed in recent numerical 
simulations [F-'7 . Moreover, enhanced spin fluctuations 
and spatial correlations in the copper oxide planes offer 
a natural framework which accounts for a depletion of 
low-energy states in the pseudogap regime of the high- Tc 
superconductors fSHTO]. As for the one-dimensional (ID) 
regime, it is dominated by spatial fluctuations [11 . The 
relevance of Umklapp scattering for the half-filled band 
leads to the absence of a bandwidth-controlled Mott tran- 
sition. However, a Mott transition can be triggered as a 
function of dimensionality. 

The aim of this Letter is to reexamine the dimensional- 
crossover-driven Mott transition in the quasi- ID Hub- 
bard model at half-filling. The subject combines many 
fascinating issues such as the breakdown of spin-charge 
separation and the binding of spinous into magnons [12 . 
These phenomena follow from a delicate interplay be- 
tween fc-space and temporal fluctuations. To capture 
the relevant physics we employ a cluster extension of 
the dynamical mean-field theory (CDMFT) [13 . In 
the CDMFT, a cluster of Nc impurities is subject to 
a dynamical effective bath simulating the effect of all 



the other sites of the lattice [14]. The ability of the 
CDMFT to reproduce the density-driven Mott transi- 
tion in the ID Hubbard model has been demonstrated 
in Refs. [15', TF. Previous CDMFT studies of the 
dimensional-crossover-driven Mott transition in weakly- 
coupled ID chains yielded ordinary open Fermi surface 
(FS) in the Hubbard model [17 and small FS pockets in 
the model of spinless fermions [18 . 

We study the Hubbard model on a strongly anisotropic 
square lattice at half- filling, 

^ = -Yl ^iAcT^j^ ^ ^ ^ ^'^^'^ -^^Ui^, (1) 

ij ,a i i,cr 

where the electron hopping tij is t {t_\_) on the intrachain 
(interchain) bonds, fi is the chemical potential and we 
set U/t = 3. In addition, we allow for a finite diagonal 
next-nearest neighbor hopping f = — tx/4 (see Fig. [T]). 
It brings about frustration in the ground state and by 
reducing nesting properties of the FS precludes long- 
range magnetic order in the T = and weak-coupling 
regime [19i[20]. Hence, finite value of t' guarantees the 
Mott transition in the thermodynamic limit. We use the 
Hirsch-Fye Quantum Monte Carlo (QMC) algorithm as a 
cluster-impurity solver and extend previous studies [17] 
to a low temperature regime. However, computational 
cost prevented us from decreasing the temperature be- 
low T = t/30 on the 8 X 2 cluster [21 . The CDMFT 
allows one to compute the single-particle spectral func- 
tion A{k,uj) = —^lmG{k,uj). Here G{k,uj) is the lat- 
tice Green's function represented in the original Brillouin 
zone. We estimate the latter by periodizing the cluster 
Green's function and applying stochastic analytical con- 
tinuation of the QMC data [22 . 

Our main results are summarized in Figs, [l] and 
|2j The control parameter t± interpolates between one 
and two dimensions and triggers the Mott transition at 
t±/t 0.18. To pin down the nature of the transition, 
continuous or first order, we plot in Fig. jl|a) = 

Efca |fr(4a^fca) ^ith Sfc = -2{t COS + ti_ COS ky) + 
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FIG. 1. (color online) (a) dF/dt± (squares) and suscepti- 
bility of interchain charge fluctuations d^F/dt]_ (scaled by a 
factor 0.1, circles); (b) spin susceptibility Xs{'i)^ and interac- 
tion vertex contribution to the d-wave pairing susceptibility 
Xp (multiplied by 100) on the 8x2 cluster. Inset shows the lat- 
tice geometry of the Hamiltonian ([T]). Parameters: [//t = 3, 
t' -t±/4 and T t/30. 
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FIG. 2. (color online) (a) Schematic evolution of the FS. The 
emergent FS forms electron pockets around k — (=b7r/2,0) 
and hole pockets around k — (=b7r/2, ±7r). The pockets 
have anisotropic distribution of the spectral weight between 
the main (solid) and ghost (dashed) vertical segments. The 
anisotropy and pockets' size grow until t^/t=0.25 when the 
ghost segments vanish and the main ones merge into the open 
FS of a quasi-lD metal (bold). The QP weight retains a 
strong Aj-dependence with a minimum at hot spots (dots), 
(b) QP weight Zfc on the main and ghost sides of the electron 
(squares) and hole (circles) pockets. 



t± cos kx cos ky. Down to the considered temperatures, 
we do not detect a jump and thus conclude that the tran- 
sition is continuous. d'^F/dt\ in Fig. [ija) corresponds 
to an interchain charge susceptibility which is greatly en- 
hanced in the vicinity of the critical coupling t^/t = 0.18. 
The open warped Fermi lines which form at sufficiently 
large t^/t > 0.24 (see Fig. [2| essentially follow from the 
topology of the tight-binding model. In the intermedi- 
ate region, 0.18 ^ t±/t ^ 0.24, we find a metallic phase 
where the FS is broken into electron and hole pockets. 
On the one hand, starting from the ID Mott insulating 
state, the occurrence of the pockets might be understood 
by taking the interchain hopping into account at the 
random-phase approximation (RPA) level [23 . In this 
context the nodal points k = (±7r/2, ±7r/2) play a spe- 
cial role since |p- = there. On the other hand, starting 
from the large t± limit, scattering off Qi = (7r,0) and/or 
q2 = (tt, tt) magnetic fluctuations could equally gap out 
the hot spots, k = (±7r/2, ±7r/2). To provide support for 
this scenario, we plot in Fig. [TJb) the cluster spin sus- 
ceptibility Xs{q) = wjo dr^ij ^''^'~'HSi{r)Sj{0)) for 
both momenta. As apparent, ID fluctuations Qi = (tt, 0) 
remain robust up to t±/t = 0.16 but are then gradually 
suppressed, giving way to dominant ^2 = (tt, tt) fluctua- 
tions, which peak at the Mott transition. Let us however 
note that in the static mean- field limit, antiferromagnetic 
order is not sufficient to reproduce the observed FS topol- 
ogy pointing towards the remnant ID Umklapp scatter- 
ing at the nodal momenta as the origin of the pockets. 

Enhanced staggered magnetic fluctuations oc Xs('7r,7r) 
give rise to pairing with a d-wave character At(r) = 

±7i[4t(^)4+5;(^) - 4ii^)4+sti^)] with the upper 
(lower) sign corresponding to S = ax icLy)^ respectively. 
The response of the system in the particle-particle chan- 



nel is best seen in the pairing interaction vertex Xp [24] . 
It is obtained from the full pairing susceptibility Xp = 
Jq^ (At(r)A^-(O)) by subtracting the uncorre- 

lated contribution. As shown in Fig. [ijb), the calcu- 
lated pairing vertex Xp is enhanced at the Mott transition 
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FIG. 3. (color online) Dimensional-crossover-driven Mott 
transition as seen in the single-particle spectral function 
A{k,(jo)\ (a,b) t± = and (c,d) t±/t 0.18. In panel (b), 
the spinon (holon) branch corresponds to the lower (higher) 
binding-energy peak, respectively. 
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FIG. 4. (color online) Low-energy part of the spectral func- 
tion A{k,oj) at t±/t — 0.2. Vertical bars track the position of 
the peaks constituting: (a) electron and (b) hole FS pockets. 



which confirms the magnetic-pairing scenario. 

The schematic evolution of the FS surface shown in 
Fig.[2|a) stems from the calculation of the single-particle 
Green's function. In Fig. |3] we show the dimensional- 
crossover- driven Mott transition as seen in the evolu- 
tion of the single-particle spectral function A{k, u). Most 
noteworthy features in the ID limit shown in Fig. [3ja,b) 
are: (i) a well defined single-particle gap at fc = (7r/2,0) 
and equivalent points; (ii) signatures of spinon and holon 
branches especially in the vicinity of fc = (0,0) [25 , and 
(iii) backfolding of the energy bands around k = (7r/2, 0) 
and the equivalent points. Concerning (ii), already when 
t^/t = 0.15 the intensity of the spinon excitation is no- 
ticeably reduced but nevertheless one can distinguish two 
peaks [26^. As shown in Fig. [3]^c,d), they are superseded 
by a single peak with a broad shoulder at the Mott tran- 
sition. Finally, a broad quasiparticle (QP) peak is re- 
solved at t^/t = 0.2 [26 . This remnant aspect of the ID 
physics is captured in approaches starting from a frac- 
tionalized spectral function in the ID limit and treating 
the interchain hopping at the RPA level [23", "27 . As 
for (iii), the pocket emerges when one of the backfolded 
bands intersects the Fermi energy. This defines a main 
and ghost side of the pocket which we can characterize 
with the magnitude of the QP residue Z^. We extract 
this quantity by fitting the data to the Lorentzian form 
^(*^'^) - ^ (c-.f+S^+r^ plot it in Fig.gb). At the 
exception of Mott transition at t^/t = 0.18 where the 
pockets shrink and become very thin, strong anisotropy 
along the pockets is evident. We illustrate this in Fig. |4] 
by showing the low-frequency part of Aik^uo) across the 
pockets at t^/t = 0.2. The two features — one with a 
large and the second one with a small QP weight cross- 
ing the Fermi level are part of the electron [Fig.|4ja)] and 
hole [Fig. [4]^b)] pocket. The broken FS is unrelated to a 
specific ladder geometry of the 8x2 cluster and it is also 
found on the square 4x4 cluster [28^. 



FIG. 5. (color online) Real part of the zero-frequency Green's 
function around momentum = (7r/2, 0) in the : (a) insulat- 
ing Mott phase and (b) metallic phase. Insets show the low- 
frequency imaginary part of the self-energy at: k — (7r/2,0) 
in panel (a) and k indicated by the capital letters in panel 
(b). G{k, 0) for t^/t = 0.2 was rescaled by 0.5. 



The reconstruction of the FS is governed by the topol- 
ogy of the real part of the zero-frequency Green's func- 
tion G(fc, 0) [29^31]. In the Fermi liquid theory, G(fc, 0) is 
positive (negative) inside (outside) the FS, respectively, 
and changes sign by going through a pole. This con- 
trasts with the Mott insulator in which G(fe, 0) changes 
sign in momentum space passing through a zero as a re- 
sult of a diverging lattice self-energy E(A:, ic<j.,n). We ex- 
tract the latter from the Dyson's equation G~^(fc, iuJrn) = 
GQ^(k^iuj^) — Yj(k^iuJm)'> where Go{k^iuJm) is the bare 
Green's function, in combination with the spectral rep- 
resentation of the lattice Green's function G{k, iujrn) = 
J duj ^^'^K The diverging behavior of E(A:, iurn) at mo- 
mentum k = (7r/2,0) in the Mott phase is shown in 
Fig. [sja). It is the coexistence of infinities and zeros 
which accounts for the anisotropy of the emergent pock- 
ets. We focus on the electron pocket and illustrate this 
in Fig. [5|b) for t±/t = 0.18. On the one hand, in close 
vicinity of the G-point with a vanishing T^{k,iuJm 0), 
two adjacent poles of G(fe, 0) result in a very thin elec- 
tron pocket. On the other hand, diverging S(fc, iuj^ 0) 
yields a zero of G(fc, 0) at the A-point. The latter reduces 
the QP weight of the nearby ghost side. The interference 
of the neighboring pole and zero becomes stronger with 
growing t_\_ and prevented us from resolving the full struc- 
ture of G(fc, 0) already dit t±/t = 0.2. Indeed, at our low- 
est temperature T = t/30 only a pole associated with the 
main side of the pocket and a broad minimum in G(fc, 0) 
is observed in Fig.jsjb). However, as depicted in Fig.[4ja) 
the ghost side remains visible in the spectral function. In 
analogy with the density-driven Mott transition in the 2D 
Hubbard model [29l [30 , we believe that the emergence 
of a large FS above t±/t = 0.24 corresponds to a simulta- 
neous annihilation of the adjacent zero and pole leaving 
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FIG. 6. (color online) Real part of the zero-frequency Green's 
function near the nodal k = (7r/2, 7r/2) point. Insets show the 
low-frequency dependence of the imaginary part of the self- 
energy at A; = (7r/2, 7r/2) (left) and the corresponding spectral 
function from bottom (t± = 0) to top {t±/t = 0.25) (right). 
For comparison, dashed line in the right inset shows A{k,(jj) 
at hp along the (0,0) — (7r,0) direction at t±/t = 0.25. 



the pole carrying a larger QP weight. 

We turn now to the nodal direction of the Brillouin 
zone. Figure [g] reveals that G(fc, 0) remains almost un- 
changed with respect to the ID regime up to the k- 
selective Mott transition at t±/t = 0.18. This agrees 
with: (i) the behavior of i;(fc, iUm ^ 0) at fc = (7r/2, 7r/2) 
which we expect to diverge in the T ^ limit up to 
t±/t — 0.2, and (ii) the gap seen in A{k,uj). On fur- 
ther increasing tj_^ G{k, 0) develops a pole-like feature. 
As a result, a broad QP peak emerges at t±/t = 0.25. 
Its weight is much smaller than that ai kp along the 
(0, 0) — (tt, 0) direction included for comparison in Fig.|6] 
The substantial variation in the scattering rate is a pre- 
cursory indication of the broken FS at smaller t±. 

Let us expand briefly on a relationship of our results 
with real materials such as quasi- ID organic (TM)2X 
salts where TM stands for TMTTF (tetramethyltetrathi- 
afulvalene) or its selena TMTSF analog [32]. Their uni- 
fied phase diagram bears similarities to that of the high- 
Tc cuprates in the sense that as a function of pressure a 
superconducting phase emerges in the proximity of the 
insulating Mott phase. The inclusion of a finite t' in 
our studies mimics the unnesting role of pressure. As 
we show, crossing over from the ID system to a higher- 
dimensional regime involves recombination of spinous 
and holons into the conventional QPs and releases the 
charge from the confinement to the ID chains. The 
pressure-induced change in the optical spectroscopy data 
on (TMTTF) 2X is interpreted as an example of such a 
charge deconfinement [33,t34i- A quantitative compari- 
son would involve studies within a model Hamiltonian 
extended by electron-phonon coupling and long-range 
Coulomb interaction to account for a wide variety of 



phase transitions observed in the (TM)2X families. Nev- 
ertheless, we find encouraging that our CDMFT studies 
within the bare Hubbard model Eq. ([T]) support accu- 
mulating experimental evidence that superconductivity 
in the (TMTSF) 2X salts is mainly mediated by magnetic 
fluctuations [35] . Finally, the signatures of the closed FS 
contours that we find could be verified in the quantum 
oscillation experiments. 

In summary, our CDMFT simulations yield a continu- 
ous quantum phase transition between a ID Mott insu- 
lating state and a 2D metallic state. On the metallic side 
the coherence temperature below which QPs form marks 
the crossover scale and vanishes at the critical point. At 
energy scales below the coherence temperature and in 
the close vicinity of the transition point, the FS topology 
shows hole and electron pockets. We attribute their ori- 
gin to the remnant ID Umklapp scattering at the nodal 
momenta. Such a mechanism can also account for the 
pockets observed in the spinless model [18 . The evolu- 
tion of the pockets with t_\_ can be understood by tracking 
the zero and poles of the single-particle Green's function. 
At energy scales above the coherence temperature, rem- 
nant features of spin-charge separation are apparent in 
the single-particle spectral function. At the two-particle 
level, the metallic state is characterized by enhanced an- 
tiferromagnetic fluctuations in the very close vicinity of 
the critical point. These magnetic fluctuations act as 
a glue for paring correlations with a d-wave character. 
On the insulating side, the crossover scale is set by the 
Mott gap. Below this energy scale we observe robust 
ID Mott physics: aspects of spin-charge separation are 
visible both in the spectral function and (tt, 0) magnetic 
fluctuations which remain intact. Further work aimed at 
investigating the finite temperature consequences of this 
quantum critical point is presently under progress. 
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SUPPLEMENTARY DATA TO THE ARTICLE: 
DIMENSIONAL-CROSSOVER-DRIVEN MOTT 
TRANSITION IN THE FRUSTRATED 
HUBBARD MODEL 

Figure [7| shows the single-particle spectral function 
A{k, Lj) on the 8x2 cluster for t±/t 0.15 and t±/t = 0.2 
and it complements Fig. [3] in our Letter. 

Auxiliary CDMFT simulations performed on the 4x4 
cluster indicate that the broken FS is unrelated to a spe- 
cific ladder geometry of the 8x2 cluster (c.f. Fig. [8| 
Hence, FS pockets appear to be a generic low-energy fea- 
ture accompanying the Mott transition in the quasi- ID 
limit. 
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FIG. 7. Dimensional-crossover-driven Mott transition as 
seen in the spectral function A(k,uj) within the CDMFT on 
the 8 x 2 cluster: (a) t±/t = 0.15 and (b) t±/t = 0.2. Param- 
eters: U/t = 3, = -t±/4 and T = t/30. 
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FIG. 8. (color online) Low-energy part of the spectral func- 
tion A{k,u) at t±/t = 0.2 within the CDMFT on the 4 x 4 
cluster. Vertical bars track the position of the peaks consti- 
tuting: (a) electron and (b) hole FS pockets. Parameters: 
U/t = 3, = -t±/4 and T t/30. 



